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We discuss the formation of antihydrogen atoms (H) in an ultra-cold positron-antiproton plasma. For 
positron densities n p > 10 8 cm -3 the characteristic formation time of stable H is determined by collisional 
relaxation of highly excited atoms produced in the process of 3-body Thompson recombination. Relying 
on the mechanisms of "replacement collisions" and " transverse collisional drift" we find a bottleneck in 
the relaxation kinetics and analyze the physical consequences of this phenomenon. 



Recently, the bound state of an antiproton with a 
positron, i.e. antihydrogen (H), has been produced at 
relativistic energies jlj. Succesful experiments on accu- 
mulating up to 10 5 antiprotons in an ion trap at tem- 
peratures T ~ H initiated the ideas to create and 
study ultra-cold H. As discussed in ]|-||], antihydrogen 
at sub-Kelvin energies is an ideal system for testing CPT- 
invariance and the weak equivalence principle by gravity 
and spectroscopic experiments. 

Cold antihydrogen is expected to be produced by merg- 
ing cold plasmas of positrons and antiprotons H. The 
attention is now mainly focused on the formation of H 
through the process of Thompson capture in 3-body col- 
lisions 
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Neutral (anti) atoms H* produced in this process have 
binding energies e ~ T and a characteristic size of the or- 
der of the Thompson radius Rt = e 2 /T (e is the positron 
charge). They subsequently undergo collisional and ra- 
diative relaxation and are supposed to acquire binding 
energies larger than the threshold energy of ionization 
by the electric field at the ends of the ion trap. For real- 
istic fields (10 — 100 V/cm) the threshold energy £o varies 
from 20 to 60K, which corresponds to the principal quan- 
tum number no ~ 100—50. We will use the term "stable" 
for atoms with binding energies e > eo, as they leave the 
ion trap without ionization and can be further captured 
in a neutral atom trap. 

We will discuss the case of an ideal plasma where T ^> 



e 2 n|!/ 3 and, hence, we have a small parameter 



n e Rir < 1. 
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Then the event rate of the Thompson process, an p n 2 (n p 
and n e are the densities of antiprotons and positrons), is 
determined by the rate constant 
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with m being the positron mass. The numerical coef- 
ficient A ranges from 0.8 in zero magnetic field (see, 
e.g., [§) to 0.08 in infinitely high fields g. The rate 
rapidly increases with decreasing T and for sufficiently 



high positron densities n ~ 10 8 cm -3 already at temper- 
atures of a few Kelvin the Thompson process becomes 
very fast. However, the formation of stable (e > eo) 
ultra-cold H is governed by the relaxation kinetics. Thus, 
irrespective of the ratio between the corresponding re- 
laxation time t and the characteristic time of Thompson 
capture tt ~ 1/cm 2 , the former is the key parameter of 
the problem. 

For high magnetic fields (B > 5 T) and low tempera- 
tures (T ~ 4 K) used in antiproton experiments |^||[|, 
the size r of atoms with binding energies e < Sq, formed 
in the Thompson process and undergoing relaxation, sat- 
isfies the inequality 

R T >r^>r B , (4) 

where tb = vt/^b is the positron Larmor radius, vt = 
\JT jm the thermal velocity, and u>b = eB/mc the Lar- 
mor frequency. Under the condition (Q) the relaxation is 
strongly influenced by the field. At experimental positron 
densities the formation of stable H by collisional relax- 
ation dominates over radiative deexcitation. Considering 
both the mechanisms of "replacement collisions" || and 
"transverse collisional drift" || we find a bottleneck in 
the relaxation kinetics and discuss the physical conse- 
quencies of this phenomenon. 

As we consider the relaxation of highly excited neu- 
tral atoms with principal quantum numbers n > no ^ lj 
the system can be treated classically. The criterion (H) 
allows us to use the drift approximation (see e.g. [ p4[ ) 
for the motion of a bound positron: In the plane trans- 
verse to the magnetic field the positron moves in a cir- 
cle of radius Tb, and the center of the circle (oscillat- 
ing along the magnetic field direction z in the Coulomb 
potential well U(p,z) — — e 2 / ^ z 2 + p 2 ) E x B drifts 
around the antiproton at a fixed transverse distance p. 
This approximation assumes that the drift frequency 
ujd ~ ce/Bp 3 <C ojb and, hence, 
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The condition (Q) ensures that rs <C Rb <C Rt- In high 
magnetic fields the distance Rb is significantly smaller 
than the characteristic size of the atom with the binding 
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energy e « £q an d, hence, the formation of stable fi can 
be well described by using the drift approximation for 
the positron motion. 

Relaxation of highly excited atoms with the transverse 
size p satisfying the condition (||) proceeds through col- 
lisions with plasma positrons and involves two leading 
mechanisms. The mechanism of "replacement collisions" , 
first found in numerical calculation ||, assumes that 
an incident plasma positron has an impact parameter 
smaller than p and becomes bound with the antiproton 
in the course of the collision, whereas the initially bound 
positron becomes unbound and carries away the increase 
of the binding energy. The latter is large (~ e 2 /p) and, 
hence, replacement collisions provide a fast relaxation of 
the positron longitudinal motion in the Coloumb well 
U(p,z). Already after a few replacement collisions the 
longitudinal motion comes to thermal equilibrium with 
the background plasma and the positron oscillates in the 
bottom of the Coulomb well, with the characteristic am- 
plitude zt ~ (p 3 /Rt) 1 / 2 -C p- Hence, the atom size 
r rs p, and the binding energy is related to the trans- 
verse size by e = e 2 / p. 

The most important consequence of replacement col- 
lisions is the increase of the binding energy of the atom 
and, accordingly, decrease of the atom size p. The fre- 
quency of replacement collisions is ~ n e p 2 VT and the 
change of the size in each collision is of order p itself. 
Therefore, the time dependence of the size can be found 
from the equation 



P 



r) 



-n e v T p 



For p, p treated as ensemble averages the validity of 
Eq.(Q) at times t ^> tr — l/n e VTR T can be proven by 
direct solution of the kinetic equation. 

The rate of replacement collisions is independent of 
the magnetic field. These collisions provide efficient re- 
laxation of excited atoms with large p (small binding en- 
ergy) even in infinitely high fields. At finite B and suffi- 
ciently small p, but still satisfying the condition (||), there 
is another relaxation mechanism, "transverse collisional 
drift" resulting from far positron-atom collisions with im- 
pact parameters larger than Rt §. Due to such colli- 
sions a bound positron undergoes diffusion in the plane 
perpendicular to the magnetic field. The corresponding 
diffusion coefficient is the same as that for free plasma 
positrons due to far positron-positron collisions: 



TR 



(7) 



where A = log (n e i?^) -1 / 3 is the Coulomb logarithm. 
The principal difference originates from the fact that the 
longitudinal motion of bound positrons is coupled to the 
transverse one, whereas for free positrons the coupling 
is absent. Therefore, collisions between free positrons 
do not change their longitudinal velocities and, hence, 
do not contribute to the mobility pT|. On the contrary, 



far positron-atom collisions change both the longitudinal 
velocity of a bound positron and its distance from the 
antiproton. This leads to a non-zero mobility coefficient 
of bound positrons, b± = D±/T [[|. Thus, there will be 
a collision-induced drift of the bound positron towards 
the antiproton, with the velocity 
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The drift increases the binding energy of excited atoms 
and, in contrast to replacement collisions, provides effi- 
cient relaxation at small p. 

On the basis of the two described mechanisms, we can 
now draw the picture of the kinetics of collisional re- 
laxation. For p ~ Rt the ratio of the rate of replace- 
ment collisions to the rate of transverse collisional drift 
is ~ (Rx/rs) 2 - Hence, highly excited atoms with binding 
energies ~ T, formed in the Thompson process ([TJ), first 
undergo relaxation due to replacement collisions. As fol- 
lows from Eqs. (^|) and (||), the rate of replacement colli- 
sions decreases with increasing binding energy e (decreas- 
ing p) and at p ~ pb = Rb(Rt/Rb) 2 ^ 5 "C Rt becomes 
equal to the rate of transverse collisional drift. Further 
increase of e occurs due to the transverse drift, and the 
relaxation rate increases with e. Thus, the total relax- 
ation rate is minimum at p « pb, and this point actually 
represents a bottleneck in the relaxation kinetics. The 
ratio pb/ Rb 3> 1 justifies the use of the drift approxima- 
tion for p ~ pb- The characteristic time of the formation 



(6) of atoms with size p > Rb is given by r = J T dp/p, 



where p 
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p( d \ The main contribution to the 



integral comes from distances p ~ pb, and r becomes 
independent of p: 
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In the range of fields and temperatures relevant for 
producing antihydrogen in ion traps ^|), the bottleneck 
binding energy Sb — e 2 / pb oc (BT) 2 / 5 is smaller or of the 
order of the threshold energy £q for typical experimental 
conditions. Therefore, the formation of stable H will be 
determined by the characteristic time r (|^). It is impor- 
tant to emphasize that possible decrease of temperature 
or change of the magnetic field in ion trap experiments 
can not appreciably reduce r. One can also check that 
the characteristic time of radiative relaxation of atoms 
with binding energies e < Eq is much larger than r. 

Neutral atoms are produced in the inner spatial region 
of the ion trap, where the electric field is scree ned by the 
plasma, and move freely with velocities v a ~ y/T/M ( M 
is the atom mass) . If the size L of the inner region is such 
that the transit time tj^Ct, then the binding energy of 
atoms which reach the ends of the trap is smaller than 
Eq. The atoms undergo ionization by the electric field 
and return to the plasma as positrons and antiprotons. 
In this case the outcome of antihydrogen from the trap 
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is negligible. Promising is the opposite case T£ ^> r, 
where the atoms reach the ends of the trap after acquiring 
binding energies e > eo and escape from the trap_without 
ionization. Then the rate of producing neutral H by the 
ion trap is determined by the rate of Thompson capture 
N p /tt (N p is the number of antiprotons). For realistic 
parameters of the ion trap experiments n e = 10 8 
cm" 3 , T = 4 K, and B = 10 T, we have r T w 1CT 4 
s, t w s, and v a ~ 2 x 10 4 cm/s. The condition 

tl/t S> 1 requires L ~ 1 cm. As the ratio tl/t cx 
T~ 4 / 5 , and t t cx T~ 9 / 2 the decrease of temperature to 1 
K allows to have L of several millimeters and reduces tt 
by more than 2 orders of magnitude. 
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